Phase diagram of the three-band half-filled Cu-O two-leg ladder 
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We determine the phase diagram of the half-filled two-leg ladder both at weak and strong coupling, 
taking into account the Cu d x 2_ y 2 and the O p x and p y orbitals. At weak coupling, renormalization 
group flows are interpreted with the use of bosonization. Two different models with and without 
outer oxygen orbitals are examined. For physical parameters, and in the absence of the outer 
oxygen orbitals, the D-Mott phase arises; a dimerized phase appears when the outer oxygen atoms 
are included. We show that the circulating current phase that preserves translational symmetry 
does not appear at weak coupling. In the opposite strong-coupling atomic limit the model is purely 
electrostatic and the ground states may be found by simple energy minimization. The phase diagram 
so obtained is compared to the weak-coupling one. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h, 74.20.Mn 
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I. INTRODUCTION 

Considerable effort has been expended to explain the 
phase diagram of the high temperature superconduc- 
tors and the associated pseudogap phenomenon within 
the framework of quantum critical and competing or- 
der pictures. Among the several candidates for or- 
ders that could compete with superconductivity are the 
"circulating current (CC) phase, "i^ the "staggered flux 



(SF) phase 



"3,4,5,6 



"d-density wave (DDW) phase 



the "spin-Peierls phase," 



» 7.8 

"stripes".^ or the "quan- 
tum liquid crystal phase, and the antiferromag- 
netic phased The CC phase, like the SF phase, breaks 
time reversal symmetry and is characterized by circulat- 
ing currents which produce local orbital magnetic mo- 
ments. An important difference between the two, how- 
ever, is that the CC phase preserves translational sym- 
metry while the SF phase breaks it. 

It is a nontrivial task to ascertain phase diagrams of 
strongly correlated systems in two dimensions due to 
the lack of accurate, systematic, nonperturbative meth- 
ods. Mean-field type approximations always favor from 
the outset particular types of orders. Exact diagonal- 
ization is constrained by small system size. Quantum 
Monte-Carlo methods have the notorious fermion sign 
problem. On the other hand, there are reliable nonper- 
turbative methods available in one dimension such as 
bosonizationi£iiLi& and the density matrix renormaliza- 
tion group (DMRG) methodi 19 i 20 i 21 By studying ladder 
systems one can take a first step toward investigating 
some of the theoretical ideas that have been conjectured 
for two-dimensional systems, while remaining in an es- 
sentially one-dimensional setting. Moreover, studies of 
ladder systems have been strongly stimulated by syn- 
thetic compounds with the ladder structure. For ex- 
ample, Sri4Cu2404i exhibits superconductivity upon Ca 
doping and the application of pressured 

Soon after the discovery of high-temperature super- 
conductivity, one-band models such as the Hubbard or 
t-J model were suggested as the simplest microscopic 



systems that could capture the interesting physics of the 
Cu-0 system^ Indeed, most work on ladders has fo- 
cused on one-band models that consider only the Cu 
d orbital, 24 : 25 ! 26 . 27 ! 28 : 29 ! 30 . 31 ! 32 ! 33 ! 34 ! 35 ! 36 However, there 
are a number of theoretical studie a 37 ! 38 that suggest that 
neither the two-dimensional Hubbard nor t-J models 
have a superconducting ground state. It is also impos- 
sible to capture some other types of potential order like 
the CC phase within a one-band model as the circu- 
lating current pattern in the CC phase involves the O 
orbitals in a crucial wayi 3 ^ A more complete model that 
incorporates these orbitals is the three band or Emery 
mofjgljijISfii The three-band model takes into account 
both strong correlations and the hybridization of Cu and 
O orbitals. 

Relatively little work has been done on three-band lad- 
der models so fariiSi 4 ^ In this paper we obtain phase 
diagrams of the half-filled Cu-0 ladder both at weak 
coupling and at strong coupling. At weak coupling, 
we examine ladder systems of two different geometries: 
with and without outer oxygen orbitals. Renormaliza- 
tion group (RG) flows are interpreted with the use of 
bosonization followed by semiclassical energy minimiza- 
tion. Several gapped phases arise depending on the par- 
ticular values of the parameters. Within the physically 
relevant region, the D-Mott phase, which upon doping 
has a tendency toward d-wave superconductivity, arises 
in the absence of the outer oxygen sites while a spin 
Peierls phase occurs when outer oxygens are added. We 
find that the CC phase does not appear at weak cou- 
pling in the Cu-0 ladder system. At strong coupling, 
the model is purely electrostatic. Treating the hopping 
term as a weak perturbation, it is then possible to make 
connections between the weak and strong coupling phase 
diagrams. 

The paper is organized as follows. In Sec. |Hj we de- 
fine the models and discuss the weak-coupling continuum 
limit. We present the RG flow equations in Sec. IIIII In 
Sec. IIVI we bosonize the Hamiltonian and the various 
local order parameters of interest. From these, in Sec.lVI 
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(a) Five orbitals in a unit cell. 
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(b) Seven orbitals in a unit cell. 



FIG. 1: Schematic diagrams of the Cu-O ladder, (a) The 
five orbital case: There are two Cu d orbitals and three O p 
orbitals in a unit cell. Hopping integrals t dp and t dp are be- 
tween Cu d and O p; t pp is the hopping between the oxygens. 
Ud and U v are on-site Coulomb interactions, and V dp , 
and Vp P are nearest-neighbor Coulomb interactions, (b) The 
seven orbital case: There are two extra outer oxygen orbitals 
in a unit cell labeled 6 and 7. 



the phase diagrams of the Cu-0 ladders are established 
at weak coupling. The volumes of the different phases 
in the high-dimensional parameter space is quantified by 
means of Monte-Carlo sampling. In Sec. IVII we present 
the phase diagram at strong coupling and compare our 
findings to those at weak coupling. Our results are sum- 
marized in concluding Sec. IVIII Some technical details 
are relegated to the Appendix. 



bital case with two Cu d x 2_ y 2 orbitals (i = 1,2), two O 
p x orbitals (i — 3, 4), and one O p y orbital (i — 5) in a 
unit cell [see Fig. D^a)] an d (2) the seven orbital case in 
which in addition to the five orbitals there are also two 
outer oxygen p y orbitals (i — 6,7) [see Fig. QJb)]. At 
half-filling, on the average, there are eight electrons in a 
unit cell in the five orbital model and twelve electrons in 
the seven orbital model. 

The Hamiltonian we investigate in this paper is given 

by 



H — H + Hi, 



(2.1) 



where 



H = 5>n;(o:)- £ [t% tra cUx)c ja (x) 



(2.2) 



is the tight-binding Hamiltonian and 

H * = E T^W M*) - 1] + E [V<r tra n t (x) nj (x) 



+V?J lter n i (x)n j (x-l)] 



(2.3) 



is the Coulomb interaction. Each site is labeled by two 
integers (i,x) where "i" labels atoms within a unit cell, 
and "x" identifies different cells. Operators c ia (x) and 
c\ a (x), respectively, annihilate and create either Cu d 
electrons (i = 1,2) or O p electrons (i — 3,4,5 in the 
five orbital case and i = 3, 4, 5, 6, 7 in the seven orbital 
case) in a:th unit cell with spin a (a =T,I)- The num- 
ber operator rii(x) = J2a c la( x ) c ia-( x ) counts electrons 
at site (i, x). 

Consider now the hopping amplitudes between neigh- 



boring Cu and O sites, t 



dp 



and tj p , 



as well as between 
nearest O and O sites, t . Since there is no particle- 
hole symmetry, the sign of each hopping term is relevant. 
Signs of the various hopping parameters arc determined 
by the symmetry of the orbitals. By choosing appro- 
priate phases for each orbital, the hopping parameters 
t -"* ra and t™ ter in Eq. I|2.2|) for the five orbital case may 
be all be taken to be positive 

-II 



t' 



dp 



= ( -dp 



if (t,i) = (l,3) or (2,4) 
if (i,j) = (1,5) or (2,5) 
if (M) = (3,5)or(4,5) 
otherwise, 



(2.4a) 



dp 
tpp 





if (i,i) = (l,3) or (2,4) 
if (i,i) = (3,5) or (4,5) 
otherwise. 



(2.4b) 



II. MODEL 



We study half-filled Cu-0 two-leg ladders with Cu 
d x 2_ y 2 and O p x and p y orbitals as shown in Fig. ^ 
We consider two different geometries: (1) the five or- 



Hopping parameters in the seven orbital case follow a 
similar pattern. The on-site energy of an electron on 
a Cu or O site is given by or e p respectively. For 
guidance we choose on-site energies and hopping ma- 
trix elements as extracted from a density functional the- 
ory (DFT) calculation for YBCO by Andersen et 



e = e c i — e p = 3.0 eV, td v — 1.6 eV, and i pp = 1.1 
eV. Although the precise values of these parameters for 
real ladder compounds will differ from those of the full 
two-dimensional problem studied by Andersen et al., for 
concreteness we use the same values here for the ladder. 



.11 



and 



Thus for example the hopping along the legs t dp 

along the rungs tj[ p are taken to be the same, though 
this is not required by symmetry. 

Turning now to the Coulomb interactions, the on-site 
energies {7, take two different values, one for Cu atoms, 
the other for the O atoms 



Ui = 



u d in = 1,2 

U p if* = 3, 4, 5. 



(2.5) 



Nearest-neighbor interactions are 



V- 



inter 



V 

pp 





v 

pp 







if (i,j) = (l,3) or (2,4) 
if (i,j) = (l,5) or (2,5) 
if (i,j) = (3,5) or (4,5) 
otherwise, 



if (z,i) = (1,3) or (2,4) 
if (i,j) = (3,5) or (4,5) 
otherwise. 



(2.6a) 



(2.6b) 



The case of seven orbitals again follows a similar pattern. 

The Hamiltonian has the usual U(l)xSU(2) global 
charge/spin symmetry. Furthermore, it is invariant un- 
der lattice translations, time reversal, parity, and chain 
interchange operations. As stated above, there is no 
particle-hole symmetry in this system in contrast to the 
one-band model. Since the symmetry differs, the RG 
equations also differ, as discussed below in Sec. IIIII 



A. Band Structure 

We first diagonalize the quadratic, noninteracting por- 
tion of the Hamiltonian Hq to obtain the band structure 
of the Cu-0 ladder. The band structure is shown in 
Fig. |21 Evident in Fig. [21 is the fact that only two bands 
cross the Fermi level so we can analyze both ladders 
(with and without extra oxygens) in a similar manner. 
Note, however, that for very different tight-binding pa- 
rameters different numbers of bands may intersect the 
Fermi level. In the two band case, it can be easily shown 
that fcpi + fci?2 = 7r at half filling, where kFn > is 
a Fermi momentum defined by e n {kFn) = fJ- and \x is 
chemical potential. However, since there is no particle- 
hole symmetry, the Fermi velocity VFn = 9e n /dk\k=k Fn 
at the two Fermi points differs even at half filling: vpi ^ 

VF2- 

The band operator ip na (x) is a linear combination of 
the c la (x) 



(2.7) 
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(a) Five orbital case 

E(eV) 















-71 \^ 








-8 





(b) Seven orbital case 

FIG. 2: Energy spectrum of the quadratic portion of the 
Hamiltonian in momentum space, (a) Five orbital case: 
Three low-lying bands are completely filled and only the up- 
per two bands intersect the Fermi level. The Fermi velocities 
at h,F\ and kF2 differ even at half filling unlike the one band 
model. We use e d — e p — 3.0 eV, t|L = t dp = 1.6 eV, and 
t pp = 1.1 eV per the DFT calculation of Andersen et al. (see 
Ref. I43I) . (b) Seven orbital case for the same parameters as 
the five orbital case. 



where n — 1, 2, 5, (6, 7) is a band index and i — 
1, 2, 5, (6, 7) is a site index in a unit cell. It is always 
possible to choose the matrix element a in to be real since 
the Hamiltonian is time-reversal invariant. Furthermore, 
due to the symmetry under chain exchange, we have the 
following relations between a in : 



(2.8) 



4 



At weak coupling it suffices to focus only on the four 
Fermi points. Therefore, we will consider only the upper 
two bands (n = 1,2) in the subsequent calculation, since 
the completely filled bands are not active. We linearize 
the bands around the Fermi points and decompose the 
band operators into right and left moving fermion fields 

(VWrMiZ ff ) : 



With this decomposition Hq reduces to 



(2.9) 



(2.10) 



B. Interaction Hamiltonian - Current Algebra 

We now discuss the interactions between the electrons 
around the four Fermi points. We introduce the usual 
current operators2&i4 



4? =: ^nPa^nPa : i J nP = 2 ' ^ 'nPa 'a^- 



M nP = ^nPa^/jlp/J, N Paf) = i>iP a i>: 



2P0> 



(2.11) 



where er denotes Pauli matrices, P denotes the chiral- 
ity (R or L), and : : means normal ordering (vac- 
uum subtraction.) Normal ordering signs are suppressed 
in the subsequent discussion for notational convenience. 
We also follow the notation of Balents et al. 26 to permit 
comparison between the RG equations for the one and 
three-band ladder systems in the Sec. IIIII 

Interaction Hamiltonian density Tij can be expressed 
into three parts in terms of their nature. Let TLj = 

H { P +Hf ) +Uf > . There are eight allowed interactions 
connecting left and right movers, with Hamiltonian den- 
sities Hj : 



The couplings in Eq. (|2.13|) just renormalize the veloci- 
ties of the charge and spin modes, and can be neglected 
in our second order calculation. Additional M n p and 
Np a f3 operators must be introduced to describe Umk- 
lapp processes 

- Uf = g lu {M\ R M 1L + m\ l M 1r ) + g 2u (Ml R M 2L 
+m\ l M 2R ) + g xu (Ml R M 2L + M 1R M\ L 
+M\ R M 1L + M 2R M\ L ) + gtm{N Rap N La0 
+ N R a pNl a(} ) + 9t u2 (N Ra pN L0a 



(2.14) 



At half-filling the three interband Umklapp terms (g X u, 
gtui, ~gtu2) are nonvanishing. The single-band Umklapp 
term, g nu , is nonzero only if fcf„ = tt/2. The detailed 
relationships between g and U, V are listed in the Ap- 
pendix. 



III. RENORMALIZATION GROUP 

We employ the RG approach combined with bosoniza- 
tion followed by semiclassical energy minimization to 
obtain the phase diagram at weak coupling. Compar- 
ison of the method to results of essentially exact DMRG 
calculations has shown it to be qualitatively reliable in 
cases where agreement has been checked; see for instance 
Ref.lM 

The current algebra as outlined in Balents et al?~ may 
be used to obtain one-loop RG flow equations for cou- 
pling constants. At half-filling there are eleven coupling 
constants, and the full set of equations is given by 



I i 2 3 ) ( : 

9lp = P \ 9tp + J^dta - 9tu\ - 9tuigtu2 - g t u2 ) 1 



.2 1 Jl r 2 
16' 



2 2 2 2 2 

9xp = ~9tp - J^9ta- - 9tu\ - 9tul9tu2 - g tu 2 - 9xu> 



~ ^ = 9lpJlR.JlL + 92pJ 2 pJ 2 L + 9xp(JlpJ 2 L 

+J2rJil) + 9i<jJ\r ■ JlL + 92oJ 2 R ■ J2L 
+9xa(JlR ■ J 2 l + J2R ' J\l) + 9t P (L R L L 
+44) + 9t<r(L R ■ L L + 4 • 4). (2.12) 

Six additional interactions arc completely chiral 
- Tif^ = Xi P {Ji R + JIl) + \ P {J 2R + J 2 l) 

+ ^xp{J\ R J 2R + JlL'hh) + ^la(JlR ' J\R 
+^1L ' + ^2a(J2R ' J2R + *^2L ' ^2L) 

+^ x <j{J\ R ■ J2R + Jil ' J 2 l)- (2-13) 



9i<t = [ 2g tp g t a - -^9ta ~ 4 5t«i ~ ^9tui9tu2 ) - ag( a , 



g 2a = a[ 2g tp g t a ~ T>9t<r ~ ^9tui ~ ^9tui9tu2 ) - Pgi a , 



-^■9tp9to — -^9t<j ~ ^9tui9tu2 



2 _ 2 

hu2 9xa 1 



9t P — go P 9tp + Y^9oagto — g X u(gtui — gtv.2), 



5 



in terms of dual Hermitian bosonic fields <p n(T and 9 na by 



9to 



90a9tp 

+^gxu{gtui + gtu2) 



9o P - -^goa - 2gxa ) gtcr 



gxu 



2.9, 



9ta gtul 



2.9 



tp 



gta gtui 



-^9x P 9 X It • 



gtui 



'{%gtp — -^gta)gxu — gtu2gxa — gtui 



%gxp 



gtui 



1 3 
- -jgxo + agip + Pg 2p + -^{agi a + (3g 2a ) 



{2gtp + gta/2)gxu — awsw — gtu2 ^gxp 

3 f 

+ -^gx<j + ugip + (3g 2p - - {agia + (3g 2a ) 

(3.1) 

where g t = gi/[ir(v F i + v F2 )], a = {v F1 + v F2 )/(2v F1 ), 
(3 = (v F i + v F2 )/(2v F2 ), g 0p = agi p + (ig 2p - 2g xp , and 
goa = agio + Pg 2 a — 2g xc7 - The dot indicates logarithmic 
derivative with respect to the length scale, i.e., cji = 
dgi/ds, where s = InZ. 

The set of RG flow equations obtained in this sys- 
tem is different from those obtained in one band case, 
due to the absence of particle-hole symmetry in Cu-0 
ladder. Equation i|3.1 I) is the most general form of RG 
flow equations for half-filled ladder systems with four 
Fermi points and with parity, time-reversal symmetry, 
chain interchange, and U(l) x SU(2) global charge/spin 
symmetry Note that Eq. (|3.fl) is invariant under band 
interchange. Equation l|3.f|l also correctly reproduce the 
RG equations in the SU(4) limit (see Ref.|jfor example) 
of v F1 = v F2 and g xa = g tv = g xu = gtu 2 = 0. Fur- 
thermore, previously derived RG equations for half-filled 
ladders^ with particle-hole symmetry can be recovered 
by setting v F \ = v F2 . 

We integrate Eq. <|3.f(l numerically starting from ini- 
tial values of gi(0) determined by the bare interactions as 
presented in the Appendix. In general, the couplings di- 
verge at some large length scale. We integrate Eq. (|3.I|I 
until the largest coupling constant equals 0.01, so that 
the one loop RG flow equations remain valid. We have 
checked that the phase diagram so obtained is not sen- 
sitive to the choice of infrared cut-off. 



4>n 



Pa 



/2ne 



[i{Pct>na + na )l 



(4.1) 



where e is a short-distance cutoff, and again P = R or 
L = ±1 is chirality. The bosonic fields satisfy the usual 
commutation relations 

[<p na (x),(p n > a >(x')] = [6 nrT (x), 6 n ' a > (x')] =0, 
{4>na{x),9 n 'a'{x')] = mS n>n i 5 a ^> (x - x') , (4.2) 

where Q(x) is the Heaviside function. The Klein factor 
K n(7 is introduced to ensure the correct anticommutation 
relations of the original fermionic fields; n na commutes 
with the bosonic fields, and satisfies 



(4.3) 



A canonical transformation separates bosons into 
charge and spin pieces 

{4>,o) np = + (<f>^U], 

{ct>,6)nc = ^=[(cf>,9U-(^6) nl }. (4.4) 



We also define (</>, 9) rv as the following: 
{<t>,6) rv = ±=[r(<t>,6) lv + {ct>,i 



(4.5) 



where r is either + or — . With these definitions, the 
nonintcracting part of the Hamiltonian density has the 
following form: 



Ho — 



vfi + v F2 



2tt 

vfi - vfi 
2ir 



Y^^M^A-u) + (d x e+„)(d x e_ v )}. 

V 

(4.6) 



The momentum-conserving part of the interaction can 
be written as the sum of two terms, 7ij — H. 1 } 1 ^ +71^^ , 
where 

\^A rv [{d x (j) rv ) 2 - (d x 9 ru ) 2 ] 



n 



(la) 



2tt 2 



(4.7) 



with 



.4 



- -^[gip + 92p + 2rg xp ], 



Ara = --[gia + g 2 cr + 2rg xa -]. 



IV. BOSONIZATION 

To elucidate nature of the phase, we use the Abelian 
bosonization technique to interpret the action semiclassi- 
cally. The fermionic field operator tp nPtT is first expressed 



And 



(16) 



1 



T [2tg ta cos 26>_ p cos 2(j> +a - {g Xa + g 2a ) 
(Zne) 

x cos 2(j)+ <T cos 2(f)- p — 2Tg xa cos 20+ CT cos 29 - a 
+ cos26L p (f'5 f + cos20_ (T + gjT cos20_ CT )], (4.8) 



G 





<t>+P 


</> + <7 


0- p 


0-* 


4>-a 


CDW 








tt/2 





* 


SP 


tt/2 





tt/2 





* 


SF 














* 


DC 


tt/2 











* 


JJ-lvlOtt 


U 


u 


U 


* 


n 
U 


D'-Mott 


tt/2 








* 





S-Mott 








tt/2 


* 





S'-Mott 


tt/2 





tt/2 


* 






TABLE I: Potential phases of the Cu-O two-leg ladder, and 
their pinned configurations (a '*' means that the field is fluc- 
tuating.) 

with gf = g t(y T ^9t P and f = k^kii^k^i- Finally, 
the bosonized form of the Umklapp interaction density 
reads 

2 

^/ 3) = — ^cos2(t> +p [2tg xu cos29- p + 2(g tu i 
(ZTre) 

+9tu2) cos 2<j) +a + g+ cos 2<j)- IJ + tg~ cos 20_ ff ] , 

(4.9) 

witn gt = {dim + gtu2) ± 4(5t«i - gtu2)- 

Note that since the Hermitian operator T obeys T 2 = 
/, r has eigenvalues T — ±1. Furthermore, since 
[H, f] = 0, H and T can be simultaneous ly diagon alizcd. 
We choose r = 1 in this paper. See Refs. I29ll30ll3l] for de- 
tails regarding to the subtleties of the Klein factors. For 
the sake of simplicity we also suppress the Klein factors 
in what follows. 

At half-filling the system is either fully gapped (by far 
the most common case) or gapless in all sectors (referred 
to as C2S2 - see 

Ref.|M) We do not find any partially 
gapped phases at half-filling. For the gapped phases the 
ground state configuration of the bosonic fields can be 
determined by minimizing the energy of the low-energy 
effective Hamiltonian at the end of the RG flow. 

Once the ground state configuration of the bosonic 
fields is determined, the bosonized order parameters may 
be examined to determine the physical nature of the 
ground states. We consider the same order parameters 
studied in previous work on one-band ladder systems: 
the order parameter of the (n, tt) charge density wave 
(CDW) phase, the (tt, tt) spin Peierls (SP) phase, the SF 
phase, the diagonal current (DC) phase, and the four 
quantum disordered phases (D-Mott, D'-Mott, S-Mott, 
and S'-Mott.) Order parameters expressed in terms of 
bosonic fields are as follows: 

Ccdw oc cos tfi +p sin cos <fi +cr cos0_ ff 

+ sin +p cos 9-p sin <p +CT sin0_ (T ,(4.1Oa) 

Osp oc cos 4>+ p cos 9- p sin sin 0_ CT 

+ sin <f) +p sin 9- p cos <fr +cr cos 9- a , (4. 10b) 



Osf oc cos <p +p cos 9^p cos 4> +a cos 9-a- 

+ sin 4>+ p sin 9- p sin <$>+ a sin^Lo., (4.10c) 

Odc of cos (j)+p sin 9- p sin <j) +a sin 9- a 

+ s'm<fr+ p cos 9- p cos 4>+ a cos9-a, (4.10d) 

CD-Mott oc e ie+p cos 6L p cos <p +a cos <j>-„ 

-ie t9+p sin0_ p sin 4> +a sin0 +CT ,(4.1Oe) 

Cs-Mott oc e ze+p cos 9-p sin <j> +a sin <p- a 

-ie l9+p sin 9- p cos <p +a cos <p- a .(4.10f ) 

We emphasize that Klein factors have been taken into 
account in the derivation of these expressions, though we 
suppr ess the K lein factors for the sake of simplicity (see 
Refs. 12913(1311 for more details.) Note that ©D-Mott is 
the order parameter for both D-Mott and D'-Mott phase 
(the same holds between S-Mott and S'-Mott.) These 
primed phases are half-cell translated states of the un- 
primed states since tfi+ p differs by 7t/2j2£L The bosonized 
form of the order parameters is essentially the same as in 
one band case except for the appearance of complicated 
band coefficients a^- . 

Table [I] lists the bosonic configurations for all the 
gapped phases. Note that since 4> rv and 9 rll are con- 
jugate to each other, they cannot be pinned simulta- 
neously. The nature of these phases was discussed in 
Refs. |2ll29||30ll3ll in one band context. Transitions be- 
tween the various phases have the same critical behav- 
ior as in the one-band case (see, for example, Fig. 5 in 
Ref.|U) 

Note also that the order parameter of the circulat- 
ing current phase that preserves translational symmetry 
does not appear in Eqs. I|4.10|l . This is because the long- 
wavelength part of the current is determined by gradi- 
ents of the boson fields, d x 9± p , but these cannot acquire 
non-zero expectation values at the semiclassical minima. 
Field 9- p locks into the particular values specified in Ta- 
ble [U and field 9 +p fluctuates with zero mean gradient. 
The bosonized form of the current along any link inside 
the unit cell, apart from the unimportant gradient terms, 
is given by 

jij(x) oc {-l) x tij(aaaj 2 - oaOji) 

x (cos 4> +p cos 9-p cos 4> +a cos 9- a 

+ sin (f) +p sin#_p sinc/(> +(T sin 6Lo-). (4.11) 

It is clear from this expression that the current pattern, 
if it exists, must be staggered along the leg direction 
such that it breaks translational symmetry. Hence, the 
circulating current phase cannot occur at weak coupling 
in half-filled Cu-0 two-leg ladder system. Furthermore, 
this is also true away from half-filling, at least as long as 
the Fermi level intersects only two bands. 
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D-Mott 


D'-Mott 


S-Mott 


S'-Mott 


Volume(%) 


40.78 


3.96 


12.7 


22.97 




SF 


SP 


CDW 


C2S2 


Volume(%) 


0.84 


5.63 


12.86 


0.19 



TABLE II: Volume of each phase over the entire parameter 
space determined by Monte-Carlo sampling in the five orbital 
case. We sample Ud, U p , Vj p , V d p, and V pp uniformly. 



V dp 




FIG. 3: Weak coupling phase diagram in the five orbital case: 
Up = 0.38(7 d and V pp = 0.1 VV 



V. WEAK COUPLING PHASE DIAGRAMS 

Equipped with the above considerations, we deduce 
the phase diagram of the three-band Cu-0 ladder system 
in the weak coupling regime. By sampling the param- 
eter space at random, we may quantify the volume of 
each phase. Table HT) lists the volume in the five orbital 
case. We find that the D-Mott, S'-Mott, S-Mott, and 
CDW phases occupy most of the parameter space. The 
SF phase also arises, though its volume is very small. 
However, the DC phase does not appear. 

It is also useful to consider some slices through param- 
eter space that may be relevant for real systems. In Fig. [3] 
we present the weak coupling phase diagram of the five 
orbital Cu-0 ladder for U p = 0.38[/d and V pp = 0.1 Vd P - 
In the absence of nearest neighbor Coulomb interactions 
the D-Mott phase appears for positive on-site Coulomb 
interactions; the S-Mott phase appears for negative on- 
site Coulomb interactions^ The CDW phase is found in 
the second quadrant bisecting the S-Mott phase. There 
is a small region of S-Mott phase located around 92° be- 
tween the D-Mott phase and the CDW phase. Note that 
direct transitions between the D-Mott and the CDW are 
not generically possible since it is necessary to unpin two 
bosonic field simultaneously: 9- p and 4>- a or 8- a (In the 



V dp 




FIG. 4: Weak coupling phase diagram in the seven orbital 
case: U p = 0.38Ud and V pp = 0.1Vd p . 



one-band model, for example, the SF phase mediates the 
transition between the D-Mott and the CDW phased). 
The SF phase arises next to the D-Mott phase in the 
fourth quadrant bisecting the D-Mott phase. A tiny por- 
tion of the D-Mott phase exists between the D'-Mott and 
the SF phases, which can be understood with the same 
argument used for the transition between the D-Mott 
and the CDW: a direct transition from the S'-Mott to the 
D-Mott or the SF phase is not possible since two bosonic 
fields must become unpinned simultaneously. Increasing 
V pp and holding other parameters fixed, the SF phase 
shrinks and then vanishes. Except for the disappearance 
of the SF phase, the phase diagram is qualitatively sim- 
ilar to Fig. |3 with slightly different phase boundaries. 

The phase diagram of the seven orbital case is pre- 
sented in Fig. The phase diagram is topologically 
rather similar to the five orbital case, however, with two 
notable differences: the SP phase now appears inside the 
D-Mott phase, and the SF phase disappears altogether. 
Both of these aspects might be related to the anisotropy 
between the effective exchange integrals parallel and per- 
pendicular to the legs. In a DMRG calculation on a 
two-leg Cu-0 ladder system by Nishimoto et al.^ it was 
found that in the absence of the outer O orbitals, the di- 
mensionless anisotropy ratio R — (Si ■ Sj) ruas / (Si ■ Sj)\ cg 
becomes significantly larger than one. Nishimoto et al. 
concluded that the inclusion of outer oxygens is crucial 
for a qualitatively correct description of two-leg Cu-0 
ladders. We can understand this behavior as follows. 
When R is smaller and the system is therefore more 
isotropic, the SP phase seems to be a more stable state 
than the D-Mott phase, as the D-Mott phase may be 
thought of as a product of rung singlets. For large values 
of the anisotropy ratio, the strong effective Heisenberg 
interaction along the rung makes rung singlets relatively 
strong, favoring the D-Mott phase. 
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FIG. 5: Strong coupling phase diagram of the five orbital 
case: Parameter values are e — 3.0 eV, U p = 0.38Ud, and 
Vpp — O.lVdp- Inset is an enlarged plot which is drawn from 
-20 to 20 eV. Phase SI through S6 are defined in the text. 



VI. STRONG COUPLING AND THE ROLE OF 
THE EXCHANGE ENERGIES 

In the strong coupling, atomic, limit of the extended 
Hubbard model, hopping matrix elements may be ig- 
nored, and the system is purely electrostatic. The num- 
ber of electrons on each orbital is restricted to be an 
integer 0, 1, or 2 by the Pauli exclusion principle. We as- 
sume periodic boundary conditions and take the number 
of sites to be large. The various charge-ordered ground 
states of the Cu-0 ladder may then be found by energy 
minimization. The lowest energy charge configurations 
at most double the size of the unit cell in size. Results 
for the five orbital case are shown in Fig. [SJ 

For positive U and V we find a uniform phase in which 
each Cu site is occupied by one electron. This phase 
corresponds to a Mott insulator. Superexchangc Hciscn- 
berg spin-spin interaction jH and J 1 - between nearest- 
neighbor copper spins are induced, respectively, along 
the legs and rungs of the ladder perturbatively at fourth 
order in the hopping4& In the five orbital case 



and in the seven orbital case 



(e + U d - U p + 7V dp - 8V pp f 

X {u~d + 2e + 2U d - U p + 8V dp - 16V PP ) ' 

(6.2) 

where the on-site energy difference between the d and 
p electrons is e = e d — e p > 0. With these effective su- 
perexchange interactions, degeneracies in the spin degree 
of freedom are lifted with antiferromagnetic interactions 
appearing when U d > U p > and V dp > V pp > 0. It is 
easy to see that the anisotropy ratio is larger than one in 
five orbital case within this simple calculation, qualita- 
tively consistent with the previous DMRG calculation^ 
This phase can then be mapped to the D-Mott phase 
found at weak coupling in a similar region of parame- 
ter space. In the seven orbital case, the exchanges are 
isotropic, weakening tendencies towards a D-Mott phase. 
It should be noted, however, that if we set V pp = in the 
five orbital case, the system becomes isotropic J" = J 1 - 
but the SP phase still does not arise. Another possible 
way to explain the SP phase is to consider frustration. 
For example, competition between spin exchange along 
the plaquette and along the diagonal can stablize the SP 
phase. Also, as pointed out by Kivelson et a/4£ other 
higher order interactions such as plaquette four-spin ring 
exchange interaction can be comparable in magnitude to 
the superexchange energy even at relatively large U, and 
other phases may arise due to the competition between 
superexchange and ring exchange interactions. 

For U < 0, the (ir, n) CDW phase arises in a similar 
region as that found in the case of weak coupling. It 
is robust against perturbation due to the hopping term. 
We also find other charge ordered phases that have no 
direct relationship to the weak coupling phase diagram. 
All of these charge ordered phases have tendencies to- 
ward stripes or charge separation. Phase 51, near the 
negative U axis with positive V, has two holes occupying 
the O p y orbital. The ladder breaks into two decoupled 
completely filled chains. In the S2 regime, two holes oc- 
cupy O p x orbitals on a leg in doubled unit cell. In the 
53 regime, pairs of two holes occupy the Cu orbital and 
the O p x orbitals on a leg in doubled unit cell. In 54 
and 55, holes gather around a Cu site. Finally, in the 
56 regime holes gather along one leg of the ladder, and 
the other leg is completely filled. 



4tj 

J 11 = 22 VII. CONCLUSION 

(e + U d - U p + W dp - W pp ) 2 

_}_ | ? | We determined the phase diagram of the half- filled 

Ud 2e + 2U d — U p + 4Vrf p — 8V PP J Cu-0 two-leg ladder both at weak and strong coupling. 

(6.1) At weak coupling, perturbative RG flows are interpreted 

4.4 with the use of bosonization. Due to absence of particle- 

— — hole symmetry the Fermi velocities of the two bands dif- 

(e + U d -Up + Wdp - 8V P p) fer to each otner and different RG flow equations gov- 

f 2 \ ern the three band model than in the case of the sim- 



U d 2e + 2U d — U p + 4V dp — 16V PP J ' pie one-band ladder. After bosonization, however, the 



Hamiltonian and the order parameters have essentially 
the same form as in the one-band case; it is the phase 
diagrams that differ qualitatively. By studying the Cu-0 
ladder with and without outer oxygen sites, four inter- 
esting conclusions may be drawn. First, in the physically 
relevant regime, the D-Mott phase arises when outer oxy- 
gen sites are absent. Second, the SP phase appears in 
the seven orbital case. Third, the SF phase only ap- 
pears when the outer oxygen sites are absent. The ap- 
pearance / disappearance of phases may be related to 
the anisotropy of the effective exchange energies in the 
system. Finally, we found that the CC phase which pre- 
serves lattice translational symmetry does not appear, 
at least at weak coupling in the Cu-0 two-leg ladder. 
Local currents or charge modulation, when they arise, 
always have a staggered pattern. 
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APPENDIX: RELATION BETWEEN THE 
CONTINUUM COUPLINGS AND THE LATTICE 
MODEL INTERACTION PARAMETERS 

In this appendix we discuss the relations between the 
continuum coupling constants A and g and the lattice 
Coulomb interactions U and V. Define w$ and w$ as 
follows: 

4? = UAj + V^ r \ (A.la) 



w$ = V™ ter , (A.lb) 

where V™ tra , and V^ ter are given in Eqs. (|2.5|l and 
(|2.6|) . Parameters and represent the intracell 
and intercell interactions, respectively. Double-counting 
is avoided by taking to not be symmetric. We in- 
troduce coefficients f(i,j, k, I, to, n,p, q) as 

f(i,j,k,l,m,n,p, q) 

= -a ik aua jm a jn cos[(k Fm + pk Fn )q]w^\ (A. 2) 

where i and j are site indices, k, I, to, and n are band 
indices, p — ±1, and q = 0, 1. 
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With these definitions, the continuum coupling con- 
stants A and g are given by 

Aip = 1,1, 1,1,-1, g), 

A2„ = ]T/(M, 2,2,2,2,-l,g), 

K P = 1,2,2,-1, q) + f(i,j,2,2, 1,1,-1, q) 

-f{i,j, l,2,2,l,-l, g ), 



AlCT — A2cr — 0, 

Ax* = -4^ /(*, £1,2,2,1,-1,0), 
~g lp = 2 /(*> J> X ' - 1 ' «) " X > X ' «)> 

hp = E 2 ^' * 2 ' 2 ' 2 ' 2 ' - 1 ' 9) - /(*' 3> 2 ' 2 ' 2 ' 2 ' «)' 

hp = ^/(W, 1,2,2,-1, 9) + 2, 2, 1,1,-1, 9) 
-f(i,j, 1,2, 2, l,l,g), 



ft* = -4j^ 1,1, 1,1, 1,9), 

ff 2CT = -4^/(i,j,2,2,2,2,l,9), 
fe CT = -4^/(i,i,l,2,2,l,l,9), 
hp = ]T 2f(i, 3, 1, 2, 1, 2, -1, q) - f(i,j, 1, 2, 1, 2, 1, 9), 
~9ta = -4 1,2, 1,2, 1,9), 

fc„ = X)2/(i,i,l,2,l,2,l,9), 
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gtul = J2 /(*» 3, 1> 1. 2 > 2 : 1> + /(*. i: 2 > 2 . 1. 1, 1. «), 

^ = -2 ]£/(»,/, 1,2, 2, 1,1,9). (A.3) 



indices in the unit cell. Equations l|A.3|) are the most 
general form of coupling constants for the case of two 
bands cutting the Fermi level and short range Coulomb 
interactions. By modifying or further general- 
izations can be studied. 



The sum in each of the above equations runs over all site 
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